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ABSTRACT 

We present spectrum synthesis fits to 808 K+A galaxies selected from the Sloan Digital 
Sky Survey (SDSS) and population synthesis of their spectral energy distributions, 
extending from the far UV (0.15 /im) to the mid IR (22 fim), based on the results of 
STARLIGHT code fitting to the SDSS spectra. Our modelling shows that K+A galaxies 
have undergone a large starburst, involving a median 50% of their present stellar 
masses, superposed over an older stellar population. The metal abundance of the 
intermediate-age stars shows that star formation did not take place in pristine gas, 
but was part of a dramatic increase in the star formation rates for originally gas- 
rich objects. We find no evidence for on-going QSO activity in the UV, which is 
well modeled by the emission of intermediate-age stars. We use K+A galaxies as local 
counterparts of high redshift objects to test for the presence of Thermally Pulsing AGB 
stars in similarly-aged populations and find no excess in the infrared due to emission 
from such stars, arguing that more distant galaxies are indeed old and massive at 
their redshift. All of our galaxies show significant excesses in the mid-IR compared to 
the light from their stars. We fit this ad hoc with a 300K blackbody. Possible sources 
include TP- AGB stars, obscured young star clusters and hidden AGNs, heating a 
significant dust component. 
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1 INTRODUCTION 

Most galaxies in the Universe fall into two broad classes 
in the colour-luminosity plane: a 'red sequence' of largely 
quiescent, typically early- type, old and metal rich sys- 
tems, or within a 'blue cloud' of star-forming, later-type, 
more metal poor and usually less massive objects. The bi- 
modality appears to be robust as a function of environ- 



incut and is observed even at higher redshift ( 


Stratcva et al. 


2001 


; Baldrv et alj|2006: Brammer et alj200£ 


; Muzzin et al. 


2012 


1H The origin of the bimodal distribution and how 



galaxies migrate between the blue cloud and red sequence 
(and vice versa) are topics of active current research, both 
observational and theoretical. It is clear that in this context 
galaxies in the 'green valley', crossing between the two main 
clumps in colour-magnitude space, may provide important 
clues as to the processes leading to the origin of the red 



1 The literature on this subject is vast and it is impossible to do 
justice to the issue within a short introduction; here we mention a 
few papers we find especially useful as starting points for readers 
interested in exploring this topic further 



sequence and blue cloud, and how galaxies evolve between 
them. 

In this framework, the class of objects known as K+A 
galaxies is particularly interesting. These appear to be galax- 
ies whose spectrum shows signs of a recent (within 0.2 - 2 
Gyr) episode of star formation, producing a strong Balmer 
series due to A stars, superposed over an old stellar pop- 
ulation (dominated by K giants) typical of spheroids and 
early-type spirals. Nevertheless, there is no strong [Oil] or 
Ha emission indic ating that star formation is currently not 
takin g place (e.g., ICouch fc Sharpies! 1 19871 : iQuintero et alj 
|2004 ). In many cases, the starburst appears to have com- 
prised j L significait fraction (20-60%) of the galaxy mass 
fe.g.. lKavirai et al.|[2007i : Choi, Goto & Yoon 2009 and see 
below). K+A's, therefore, are observed in a short-lived phase 
as they transit from the blue cloud to the red sequence and 
thus provide the opportunity to obtain unique insights into 
the processes leading to the morphological and spectropho- 
tometric transformation of galaxies. 

Major starbursts of the importance observed in these 
objects are often encountered in mergers and interactions 
(e.g., the Antennae) and there is in fact some morphological 
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evidence that a substantial fraction of K+A galaxies may be 
in interacting sys tems (see e.g. Yamauchi, Yagi & Goto 2008; 
lYang et al.ll2008l and references therein). While we therefore 
have some hints as to what triggered the massive starbursts, 
the real mystery is what made star formation stop so sud- 
denly that today we basically observe no traces of on-going 
star formation. 

However, a significant fraction of the current star forma- 
tion in distant clusters (where K+A g alaxies were originally 
discovered by iDressler fc Gunij[l983 ). is obscured by dust 
(Due ct al. 2002, Saintonge, Tran & Holden 2008, Dressier et 
al. 2009) suggesting that dust obscuration may explain the 
lack of observed star formation in K+A galaxies. Shioya, 
Bekki & Couch (2001) suggest that the progenitors of K+A 
galaxies may have been dusty starbursts of the e(a) type. 
It is also possible that star-forming is continuing: A-type 
stars could migrate away from their parental dust clouds 
during their lifetime, whi le more massive stars would still 



remain heavily obscured JPoggianti Sz Wu I7UUU I. tnus pro- 
ducing an A-rich spectrum. ISmail et al.l ( 19991 ') found evi 



2000), thus pro- 



dence of on-going star forma tion in their sample of 5 K+A's 
while iMiller fc Owed (|200lh detected 1.4 Ghz emission in 
2 out of 15 K+A galaxies (this may be due to an AGN or 
to synchrotron radiation fro m supernov a remnants in an ob- 
scured star forming region). iGotol |2004l set a limit of 15A/q 
yr _1 to the current star formation in a few K+A galaxies ob- 
served with the VLA. At least some K+A galaxies are seen 
to have substantial gas reser voirs comparable to those of spi- 
rals o f similar luminosities (jChang et al.l[200ll : iBuvle et alJ 
2006), suggesting that these galaxies may potentially be able 
to reinitiate star formation in the future. 

In our previous work we have found that a stacked 
FIRST (Becker, White & Helfand 1995) image of 811 
K+A galaxies selected from the SDSS Data Release 7 
l|Abazaiian et al]|2009l ; lGotoll2007l ) shows only weak 1.4Ghz 
continuum signal, most of which appears to come from 
a subsample of about 5% of the galaxies that have (5a) 
unambiguous detections in the FIRST survey (above 0.8 
mjy), approximately equally split between AGNs and star- 
forming galaxies. The median stack of the remainder of 
the g alaxies is consisten t with a purely radio-quiet popula- 
tion (|Nielsen et al.ll2012f ). This may suggest that most K+A 
galaxies are observed in the 'off' phase of an activity cy- 
cle, and that such galaxies might actually move back and 
forth between the blue cloud and red sequence until their 
gas reservoirs are finally exhausted. Such a cycle of activity, 
m odulated by mergers and 'feedback', has been proposed 
bv lHopkins et all ([2008 ) to account for the formation of a 
tight red sequence within a hierarchical merger scenario. It 
is therefore possible that K+A galaxies may be 'playing pos- 
sum' rather than having permanently migrated towards the 
red sequence following a (merger-induced?) starburst. 

K+A galaxies may also provide clues as to the nature of 
high redshift objects and their star formation history. There 
is a significant population of 1 < z < 2 galaxies that ap- 
pear to be red, quiescent, and massive. Such systems are 
difficult to explain within standard ACDM formation sce- 
narios, as they appear to b e already well-fo r med despite the 
Univ er se's relative you th |Yan et al.l |2004| ; ICirasuolo et al.l 
|2010T ). iMarastonl (120051 ) proposed that these objects 'fake' 
their old ages because they are dominated by ~ 1 Gyr old 
stars, which produce Thermally Pulsing AGB stars (TP- 



AGB), whose high luminosities and red colours mask the 
underlying blue light and th erefore yield apparently larger 
stellar masses and older ages dTonini et al.ll2009l .l2010), oven 
for young and relatively low mass objects. Present-day K+A 
galaxies of course also host older stellar populations, but 
they have experienced a major star-formation episode, in- 
volving a large fraction of their stellar mass, and may there- 
fore be adequate l ocal counterparts t o test this hypothesis 
i|Kriek et alj|20ld : Izibetti et al'll2012l ), whereas detection of 
the spectral signatures of TP- AGB stars would be very dif- 
ficult at high redshifts. 

In this paper we present a statistical stud y of the prop- 
erties of the 811 K+A galaxies studied by iNielsen et al.l 
for which we have combined SEDs obtained us- 
ing archival data from GALEX, SDSS, 2MASS and WISE, 
with the population s ynthe sis study of the SDSS spectra by 
ICid Fernandes et al.l l|2005l ). We use the combined knowl- 
edge of stellar populations, models, and spectral energy 
distributions covering the entire thermal spectrum over 2 
decades in wavelength (from 0.15 to 22 (im) to discuss the 
internal constituents of K+A galaxies, test current star for- 
mation and the nature of the starburst, and compare with 
high redshift galaxies. In the next section we describe the 
dataset, photometry and spectroscopy, and population syn- 
thesis models. We then discuss and interpret the results and 
present our conclusions and suggestions for future work. Un- 
less otherwise stated and where applicable, we use the 'con- 
cordance' cosmological par ameters deriv e d fro m the latest 
analysis of WMAP data in I Jarosik et al.1 i|201ll ). 



2 DESCRIPTION OF THE DATASET AND 
BROAD PROPERTIES 

We used a sample of 8 11 K+A galaxies s elected from the 
SDSS as described by l|Nielsen et all 12012). Th is is an ex- 
pansion of the K+A sample used by Goto! ( 2007 ) and drawn 



from the 7 th data release of the SDSS (Abazaiia n et all 



2009). Briefly, only objects classified as galaxies with a 
spectroscopic S/N > 10 per pixel are considered. The se- 
lection criteria of K+A galaxies are equivalent widths of 
H Q > -3.0A; Rs > 5.0A; and [OII]> -2.5A where emission 
lines are negative. Galaxies with redshift 0.35 < z < 0.37 
have been excluded due to contamination by the 5577A at- 
mospheric feature. 

Figure [T] summarizes the principal characteristics of the 
sample. Galaxies span the redshifts < z < 0.4 with a me- 
dian of z = 0.138 and have typical luminosities of ~ 1O 1O L0 
and total stellar masses (as given by STARLIGHT; see below) 
of ~ 3 x lO lo A/0 comparable to approximately L* galaxies, 
with a longer tail towards lower mass objects. Most objects 
are resolved and have apparent radii of a few arcsec, typical 
of galaxies of their luminosity at the appropria t e red shift. 

Of the 811 galaxies in the INielsen et ail (|2012h sa.m- 
ple, we have retrieve d spectra from the SDSS database 
l| Abazaiian et al]|2009l ). as well as the associated SDSS pho- 
tometry in ugriz. For these galaxies we then proceeded 
to obtain multi-wavelength pho tometry in the far ult ravi- 
olet from the GALEX archive ilMorrissev et al.ll2007l). in - 
frared (JHK) data from the 2MASS l|Skrutskie et al.ll2006l ) 
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Figure 1. Properties of the K+A sample used in this paper. 
From top to bottom we show the redshift distribution, distribu- 
tion of stellar masses from STARLIGHT (see below) in solar units, 
the distribution of r band luminosities in solar units, and the an- 
gular size distribution in seconds of arc as characterized by the 
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and the latest release of the UKIDSS surve}0, and finally 
fluxes at 3.4, 4.6, 12 and 22 /im from the WISE dataset 
( Wrig ht" et al.l l2010l ). We also used archival data from the 
Spitzer Space Telescope for a few galaxies observed as part 
of other programs (e.g., Lockman Hole, ELAIS). 

Given the fact that this study is entirely based on exist- 
ing archival data, from the outset we were concerned about 
possible systematic effects due to aperture matching in the 
photometry. Thus, we carefully inspected the relevant data 
bases for the measurements that would minimize any sys- 
tematic difference for the different catalogues. 

Our base photometry is the optical u, g, r, i, z from the 
SDSS, which is measured in Petrosian apertures that are 
expected to approximate the total flux of galaxies. GALEX 
photometry is calculated in similar Kron-style apertures, but 
the resolution (1.5") and point spread function (4.5" in the 
near UV and 6.0" in the far UV) are such that, given the 
sizes of most our objects (Figure [TJ, they are basically star- 
like as seen by this telescope. 

The 2MASS survey provides isophotal and extrapo- 
lated magnitudes for objects in the Extended Source Catalog 
(XSC), comprising 291 of our 811 K+A galaxies. Compari- 
son with UKIDSS Petrosian magnitudes (for 71 of these ob- 
jects which are in common with the XSC) showed that the 
isophotal estimates are more stable and closer to a 'SDSS- 
style' total magnitude, and these were adopted. Almost all 
of our galaxies are also in the 2MASS Point Source Cata- 
log (PSC). For the remaining 620 objects (in the PSC) we 
also carried out a comparison with the UKIDSS photome- 
try in the appropriate bands and eventually concluded that 
the 2MASS 4" aperture magnitude yields a closer estimate 
to the total magnitude and this was therefore employed in 
our analysis (note that this is comfortably larger than the 
Petrosian sizes of most of our sample as shown in Figure [T] 
above) . 

WISE gives magnitudes within the XSC isophotal aper- 
tures if objects are positionally matched to 2MASS and 
'adaptive' apertures for all other galaxies. We accepted the 
WISE estimate of the best magnitudes for our galaxies, as 
provided by their database. As it turned out, 12 of our 
objects have been observed by IRAC and MIPS24 with 
SPITZER. The corresponding SPITZER Kron magnitudes 
compare very well with the WISE values in the 3.4 — 3.6/ira 
and 4.5 — 4.6/xm bands, which gives us confidence that the 
WISE magnitudes are measured through apertures consis- 
tent with the rest of the bands. 

Table Q] shows the first 10 objects and their photom- 
etry; the full table is available on the MNRAS website in 
machine-readable format. Only for 98 galaxies in our sam- 
ple we have full FUV(0.15^m) to MIR (22^m) SEDs; while 
for essentially the full sample we have u' to 12/x data. 

The table shows the SDSS name (in the form hh- 
mmss.sss and ddmmss.ss) and the fluxes in GALEX (FUV 
and NUV), SDSS (ugriz), 2MASS (JHK) or UKIDSS 
(YJHK), and the four WISE bands (in this order) in units 

2 The UKIDSS project is defined in Lawrence et al. (2007). 
UKIDSS uses the UKIRT Wide Field Camera WFCAM (Casali 
et al. 2007). The photometric system is described in Hewett et al. 
(2006), and the calibration is described in Hodgkin et al. (2008). 
The pipeline processing and science archive are described in and 
Hambly et al (2008) 
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Table 1. Spectral Energy Distributions (xlO 17 ergs 1 cm 2 A x ) 







GALEX 

FUV NUV 




g' 


SDSS 








Y 


2 MASS o 
J 


- UKIDSS 

H 


Ka 


3.4/* 


WISE 

A.Gfj. 12fi 


22/j 


1 


001027.392-104341.71 


0.00 


14.57 


21.91 


46.73 


49.51 


41 


73 


35.40 


0.00 


21.21 


10.78 


7.50 


1.83 


0.93 


0.55 


0.39 


2 


001551.253-102317.92 


17.01 


11.74 


13.21 


35.22 


38.24 


31 


69 


26.42 


0.00 


16.41 


10.58 


5.92 


1.14 


0.49 


0.07 


0.00 


3 


010823.370-094356.22 


0.00 


0.00 


8.33 


17.70 


24.44 


22 


44 


19.46 


0.00 


12.44 


11.31 


5.15 


1.13 


0.44 


0.00 


0.19 


4 


011942.227+010751.63 


16.51 


24.61 


25.91 


56.62 


57.04 


48 


34 


42.05 


0.00 


32.57 


24.67 


11.38 


2.39 


0.81 


0.08 


0.00 


5 


012015.683-095920.12 


0.00 


10.17 


12.92 


34.19 


37.75 


32 


13 


27.35 


0.00 


22.10 


12.20 


3.63 


1.62 


0.60 


0.07 


0.24 


6 


014327.340-011342.55 


3.69 


6.66 


7.75 


20.41 


28.80 


26 


64 


22.75 


20.55 


16.15 


11.84 


7.26 


1.47 


0.64 


0.35 


0.37 


7 


014838.633-093259.77 


3.90 


0.00 


9.12 


24.25 


32.38 


28 


6 1 


24.63 


0.00 


16.19 


10.64 


6.16 


1.29 


0.49 


0.07 


0.00 


8 


015107.018-005636.72 


1.49 


4.93 


9.69 


24.30 


31.78 


28 


20 


24.74 


21.92 


18.45 


12.27 


6.95 


1.59 


0.65 


0.11 


0.00 


9 


020505.984-004345.25 


23.26 


29.85 


27.71 


47.72 


46.78 


39 


83 


31.64 


26.33 


19.83 


9.36 


5.98 


1.18 


0.44 


0.15 


0.00 


10 


021007.624-095431.19 


227.70 


181.10 


79.85 


126.30 


105.20 


85.36 


62.76 


0.00 


0.00 


0.00 


0.00 


1.19 


0.39 


0.17 


0.11 



Upper limits are shown as 
zeros in the table. For reasons of clarity and economy of 
space, Table [T] does not give the measurement errors, which 
are publicly available from the source web-sites. It should 
be reiterated here that the magnitudes come from a vari- 
ety of methods on data and instruments of varying reso- 
lution and therefore we assume that these individual mea- 
surements adequately approximate the total magnitude of 
galaxies without distorting their spectral energy distribu- 
tions (this amounts to saying that the apertures used are 
reasonably similar and colour gradients weak). The tabu- 
lated fluxes have been corrected for extinction as described 
below. 

We will first analyze the stellar populations of these 
galaxies based on modelling the SDSS spectra and line ra- 
tios. This information will then be used to inform our un- 
derstanding of their spectral energy distributions and the 
nature of the astrophysical sources contributing to the light 
in different wavelength regimes. 



2.1 Stellar Populations from Spectrum Synthesis 
Models 



Table 2. Age and Metallicity bins used to display the results of 
the SSP model fits. 



IStasinska et all l|200ct ) have used STARLIGHT to fit sim- 
ple stellar population (SSP) models to hundreds of thou- 
sands of galaxy spectra from the SDSS/DR7 data re- 
lea se. Full details about the m odelling process are given 
by ICid Fernandes et al.l |2005h . Briefly, the spectral fits 
are based on the Bruzual & Chariot (2003; BC03) mod- 
els with the STELIB library jLe Borgne et all I200E ), and 
the Pa dova 1994 evolu tionary tracks (|Bertelli et al.l [19941 ) 
with a I Chabrierl |2003T l IMF. The data, the model fits for 
the galaxies in our sample, and the code itself are publicly 
available from the STARLIGHT website[f| and many details 
about the use of these models for studies of galaxies with 
very weak emission lines are given in the series of papers 
by the SEAGal collaboration cited bv lCid Fernandes et al.l 
|201ll ). 

The SDSS STARLIGHT database includes 808 of 
our 811 galaxies. Therefore, the three missing ob- 
jects (J103833.180+034930.29; J143852.706+613020.53; and 
J144332.436+304807.40), for which we do have the SEDs 
and spectra, will not be included in our study. The SSP 
models fit the stellar populations using 25 different ages and 
6 metallicities. Examples of STARLIGHT fits for two 'typical' 
K+A galaxies are shown in Figure [2] that presents the ob- 



http://www.starlight.ufsc.br 



Bins 


STARLIGHT 


MARASTON- 


-mix 




BC03+Chabrier 


rhb+Kroupa 




Age t [Gyr] 




young 


t < 0.1 


t < 0.1 




intermediate 


0.25 s; t ^ 4.25 


0.3 ^ t ^ 


4.25 


old 


t > 4.25 


t > 4.25 






Metallicity Z 




metal poor 


Z sg 0.008 


Z sc 0.01 




solar 


Z = 0.02 


Z = 0.02 




metal rich 


Z = 0.05 


Z =(0.04- 


f0.07)/2 



served and fitted spectra, the residuals, and the star forma- 
tion histories resulting from the fit (ages, metallicities). 

For a variety of reasons, obtaining reliable estimates of 
the statistical uncertainties in the stellar populations de- 
rive d from ssp model fits is a notoriously diffic ult task (see 
e.g. iGonzalez Delgado fc Cid Fernandes! (|2010l ) for a recent 
discussion) and in particular STARLIGHT do es not provide 
statis tical errors. In a previous investigation (1 Melnick et al.l 
I2012T) we found that using wider age and metallicity bins in- 
creases the robustness of the models against statistical and 
systematic (i.e. variations in the input parameters) fluctu- 
ations. Therefore, we binned the populations into 3 broad 
metallicity (Z) bins and 3 age groups (t) as indicated in Ta- 
ble [2] The binned populations of the first 10 galaxies are 
listed in Table (the full table is available on the MNRAS 
website in machine-readable format). 

Figure [3] summarizes the star formation histories of 
K+A galaxies as derived from the STARLIGHT fits: 

• On average, the intermediate age stars (formed in the 
recent star formation episode and contributing to the A 
spectrum) account for close to 50% of the total stellar mass 
in K+A galaxies. This strongly suggests that the starburst 
has involve d a very significant p ortion of the initial mass of 



nas involved a very significant p ortion ol trie init ial mass ot 
the galaxy (|Kavirai et al.1 120071 : Ichoi et al.ll2009h. Analysis 
of colour gradients for K+A galaxies by Yang et al.l (|2008l ) 
and IFU spectroscopy bv lPracv et al.l (|2012l ) argues that the 
starburst tended to take place preferentially in the central 
regions of galaxies. 

• In agreement with the lack of strong [OIII] and Ha 
emission, the population of young (< 0.1 Gyr) stars is neg- 
ligible; the mean of the sample is ~0.3% and only three 
galaxies have more than 10% young stars. 

• The intermediate-age (IA) populations tend to be pre- 
dominantly metal poor or metal-rich; the median galaxy has 
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Figure 2. Left. Example of STARLIGHT modelling for a galaxy where both the luminosity and the stellar mass are dominated by the 
A-component. The larger panels show the observed (black) and fitted (red) spectra and the residuals, which clearly show the emission 
lines (that have been masked from the population fit so they appear in purple). The smaller panels show the fraction of the luminosity 
as a function of age for 3 different metallicities Z (bottom), and the fraction of the stellar mass plotted in the same way on the top 
panel. The colour-coding of metallicities (blue: metal poor — Z ^ 0.5Zq; green: solar; red: metal rich - Z = 2.5Zq) is indicated in the 
lower panel, while a number of parameters relevant for the fits are displayed in the upper panel. Notice that the EW(Ha) in absorption 
is larger than that in emission, illustrating the importance of subtracting the models in order to derive meaningful emission-line fluxes. 
Right. Similar to the left panel, but for a galaxy with a predominant old population by mass. In this example the emission lines are also 
heavily contaminated by underlying absorption. 



Table 3. STARLIGHT Populations for K+A galaxies 





SDSS 






Mass 


Av [ 


mag] 


[Nil] 
Ha 


Equi 


/alent w 


idths [A] 


Young 


Interi 




ago [%] 




Old [ %] 












10 9 M Q 


Int. 


MW 


Ha 


[out] 


[Oil] 


[%] 


poor 


solar 


rich 


poor 


solar 


rich 


1 


001027.392-104341.71 





1329 


20.9 


0.32 


0.08 


1.2 


4.9 


7.6 


2.3 


0.0 


23.4 


0.5 


37.5 


8.4 


0.0 


30.2 


2 


001551.253-102317.92 





1982 


26.9 


0.38 


0.10 


1.0 


1.1 


-0.2 


1.4 


0.2 


3.7 


0.0 


85.4 


0.1 


10.2 


0.0 


3 


010823.370-094356.22 





1834 


40.0 


0.07 


0.17 


1.4 


0.7 


0.4 


0.4 


0.1 


5.0 


0.3 


30.2 


10.1 


46.5 


7.7 


4 


011942.2274-010751.63 





0900 


24.9 


0.05 


0.05 


0.3 


1.8 


-0.2 


1.9 


0.2 


7.7 


0.9 


6.1 


3.7 


28.8 


52.0 


5 


012015.683-095920.12 





1369 


25.5 


0.03 


0.14 


2.5 


0.6 


0.2 


1.5 


0.4 


12.9 


0.9 


37.8 


11.6 


35.5 


0.0 


6 


014327.340-011342.55 





1963 


59.7 


0.20 


0.04 


1.5 


4.2 


5.2 


2.1 


0.1 


5.5 


0.4 


9.7 


26.9 


28.6 


27.8 


7 


014838.633-093259.77 





2170 


29.1 


0.16 


0.05 


1.8 


0.8 


0.3 


3.1 


0.4 


21.0 


0.0 


73.1 


0.0 


2.6 


0.0 


8 


015107.018-005636.72 





1981 


29.7 


0.10 


0.04 


1.3 


0.9 


-0.2 


0.5 


0.0 


18.3 


1.2 


68.9 


11.6 


0.0 


0.0 


9 


020505.984-004345.25 





1134 


16.8 


0.56 


0.01 


0.6 


4.5 


0.9 


3.5 


0.2 


8.4 


0.3 


25.6 


34.6 


20.8 


9.4 


10 


021007.624-095431.19 





0412 


1.9 


0.3 1 


0.01 


0.4 


4.2 


0.5 


4.6 


0.1 


12.7 


0.0 


3.6 


81.6 


2.1 


0.0 



close to 80% metal-rich IA stars, while less than 10% of the 
galaxies have more than 20% IA stars of solar metallicity. 
This supports our previous point that a rapid starburst took 
place, involving a large fraction of the pre-existing gas. Mi- 
nor mergers with gas-rich subsystems are unlikely to pro- 
vide sufficient material to account for the large fraction of 
younger stars, and would presumably dilute the existing ma- 
terial with more pristine gas. 

• We suggest that K+A galaxies must have therefore ex- 
perienced a dramatic increas e in their star form ation rates 
followed by rapid quenching (|Brown et al.|[2009h . given the 
lack of current star formation. The abundance pattern ob- 
served supports this hypothesis, as we would expect a gentler 
and more continuous distribution if star formation had been 
slowly throttled by mechanisms such as 'strangulation' or 
'harassment'. 



• About 15% of the galaxies do not contain an old stellar 
population. Otherwise, the old stellar population is predom- 
inantly metal poor. Close to 50% of the galaxies do not con- 
tain old, metal rich stars, while the median fraction of solar 
abundance stars is 20%. The original galaxies must there- 
fore have app eared as late(r) - type s pirals of apparently lower 
stellar mass. iMendel et alj (J2012) suggest that the K+A 
phenomenon is related to the growth of the galactic bulge, 
and therefore morphological transformation. 

2.2 Emission-line Diagnostic Diagrams 

Almost all the K+A galaxies in our sample show weak but 
clearly visible emission-lines of [NII]6548,6584 and Ha, and 
many galaxies also show [OII]3727,3729 and [OIII]5007. It 
is possible, therefore, to use emission-line diagnostic dia- 
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Figure 3. Stellar populations by mass in K+A galaxies. Top. 
The histograms show the distribution of young, intermediate age, 
and old stars in K+A galaxies. The sample contains only one 
galaxy with a significant young stellar population; intermediate- 
age and old stars are the dominant populations Bottom. Cu- 
mulative fractions of old and intermediate-age stars in the three 
metallicity bins discussed in the text. Metal rich stars dominate 
the intermediate-age populations, while metal poor stars are the 
most common among the old stars. The metallicity distribution of 
the intermediate-age population is heavily skewed towards metal- 
rich and metal-poor stars, while solar abundance stars are rela- 
tively rare. 



Figure 4. WHAN diagram for 787 K+A galaxies with S/N > 5 
around 660nm. The traditional boundaries separating LINERs, 
SEYFERTs, and Star-Forming (SF) galaxies are shown by the 
dotted lines (see text for more details). 



grams to gain further insight into the stellar populations 
by exploring the photoionizati on properties of the sam- 
ple. The SEAGal collaboration ijCid Fernandes et alJ[201ll : 
henceforth CFS11 and references therein) have extensively 
explored this issue on a sample of 700 thousand galaxies 
from the SDSS with weak emission lines. In particular, they 
noted that the traditional 'BPT' diagnostic diagrams that 
make use of the [Oil] , [OIII] , and H/3 lines (among others) 
are not suitable for galaxies with weak emission lines because 
these lines are usually very weak and, in the case of H/3, are 
heavily contaminated by underlying absorption. This is par- 
ticularly important for K+A galaxies where the emission 
lines are, by definition, weak while Balmer absorption lines 
are strong. 

They therefore introduced a new diagnostic diagram, 
which they call 'WHAN', that only uses [Nil] 6584 and Ha, 
whi ch are usually the strongest l ines in these types of galax- 
ies. ICid Fernandes et all (|2010l ') and CFS11 have in fact 
shown that their EW(Ha) vs. [NII]/Ha 'WHAN' diagram 
retains the diagnostic power of the more traditional BPT 
plots, but also allow to separate 'true' from 'fake' AGNs as 
we shall see below, while being less sensitive to reddening. 

Figure 0] shows the WHAN diagram for the 787 galaxies 
in our sample that have S/N > 5 in the continuum around 
[Nil] & Ha (as indicated in the top-label of the figure). 
The dashed lines separate the regions populated by LIN- 
E RS, SEYFERTS , and s tar-forming galaxies (SF) as defined 
by iKewlev et al.l J2006 ) trans ported to the WHAN plane 
bv lCid Fernandes et al |2010l ). The plot shows that, as ex- 
pected, the overwhelming majority of our K+A galaxies fall 
in the LINER domain. 
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2.3 Spectral Energy Distributions (SED) from 
0.15 to 22 

We can now analyze the full UV-optical-infrared spec- 
tral energy distributions of K+A galaxies in the light 
of the information on their stellar populations provided 
by the analysis of their spectra and the WHAN dia- 
gram. Fluxes retrieved from the archives were corrected 
for foreground Galactic extin ction using the Ay values of 
ISchlaflv fc Finkelbeinerl l|201ll ) (downloaded for each source 
using NED) and for internal extinction using the Ay given 
by the STARLIGHT models. In both cases we used the COM 
l|Cardelli. Clayton, fc Mathvslll989l ) extinction law for the 
Milky Way in view of the high metallicity of our galaxies. 
Figure [5] presents sample SEDs for two of the K+A galaxies 
with full FUV to MIR data. 

The STARLIGHT models use BC03 but one of our aims in 
this study is to verify whether TP-AGB stars are needed to 
reproduce the spectral energy distributions of 'young' galax- 
ies (whether K+As or objects in the distant universe) and 
we therefore test our observations against the SSP models of 
Maraston (2005, hereafter M05), which explicitly treat this 
stage of stellar evolution (see Introduction). We used the 
M05 SSP models for the Kroupa IMF, which provide spec- 
tral energy distributions for 6 different metallicities ranging 
from ~l/180th to 3.5 times solar (-2.25 to +0.67 dexfl and 
two different horizontal-branch (HB) morphologies, chosen 
heuristically to account for the second parameter problem 
(bhb for blue and rhb for red horizontal branches). 

To simplify the analysis, and because at this stage we 
are mostly concerned with the ensemble properties of K+A 
galaxies rather than the detailed modelling of individual ob- 
jects, we averaged together the models into three age bins 
and three metallicity bins matching as close as possible the 
ranges spanned by STARLIGHT as shown in Table [2] 

We then used the stellar mass fractions for each 
STARLIGHT bin to generate our model SEDs as: 

£-=i Si 

where for each galaxy Si are the STARLIGHT mass fractions 
for the six age and metallicity bins described in Table [2] (in- 
cluding only the intermediate-age and old populations since 
the fraction of young stars is negligible); Mi (A) are the cor- 
responding MARASTON bins defined in Table [2] 

Notice that since we are not including the young pop- 
ulations, "}2 Si ^ 100% although for 95% of the galaxies the 
sum is > 96% and only for 4 galaxies is the sum < 80%. 
Throughout the paper we will refer to the resulting DJl\ as 
the MARASTON-mix for each galaxy, for which the constant 
A is fixed matching the model to the observed fluxes in the 
SDSS r'-band (0.626^) at the corresponding redshifts. 

Figure [5] shows examples of a very good fit (left) and 
a rather poor fit (right). Since the differences in HB mor- 
phology are only relevant for old stars, we chose to use the 



4 We did not use the [M/H]=— 2.25 models, however, in order 
to maximize overlap with the STARLIGHT models that stop at 
[M/H]=-1.35. 



rhb models that cover the full range of agesQ The legends 
also indicate the fraction of the total stellar mass that was 
included in the composite model spectra for each galaxy. 

Figure [5] (see also Figs. [TTJ] and I12p clearly shows that 
our Maraston-mix SEDs reproduce quite well the observed 
SEDs from the FUV to the optical (including z'), but deviate 
systematically towards the red; there is a clear excess in the 
IR, which we ascribe (see later) to the treatment of TP-AGB 
stars, and a significant excess over the expected blackbody 
behavior at longer wavelengths, where a hot dust component 
appears to be needed to account for the observations. 

Determining the temperature of this hot dust compo- 
nent on the basis of only the WISE data, however, is rather 
uncertain so we searched the IRAS and SPITZER data bases 
to extend the wavelengths coverage in the mid- to far-IR, 
and found SPITZER data for 12 sources in the 3.6^ - 24^ 
range, although only a few of these galaxies have complete 
wavelength coverage in this range. 

The data for one such galaxy are plotted in the left 
panel of Figure[S]as blue triangles. Using the SPITZER data 
we found that Tdust ~ 300°K provides a good fit to the 
observations of these 12 galaxies, and we therefore adopted 
the same temperature for the hot dust component of all 
the galaxies in our sample, although for individual cases, 
such as the one shown in the figure, a different temperature 
(in this case ~400K) would have provided a better fit. The 
uncertainties in the Maraston models discussed below, and 
the fact that only slightly more that 25% of our galaxies are 
detected at 12 and 22 microns, however, preclude leaving 
Tdust as a free parameter. 



3 DISCUSSION 

Analysis of the spectra of K+A galaxies shows that these ob- 
jects appear to consist of originally late-type galaxies that 
have experienced a starburst involving a large fraction of 
their stellar mass. We have modeled their stellar popula- 
tions on the basis of these spectral fits and we can now 
compare with the observed SEDs to test models for the ori- 
gin of the UV flux, the presence of TP-AGB stars, and the 
contribution from hot dust. 

We assume that the stellar populations deduced from 
the SDSS spectral fibers can be extrapolated to the galax- 
ies as a whole. Only integral field spectroscopy can resolve 
this difficulty, but it is not generally available for our ob- 
jects. It seems reasonable to assume that the SDSS spec- 
tral fibers are well centered on the brightest parts of our 
relatively compact galaxies, but we can not exclude the pos - 
sibility of radial population gradients (|Pracv et alj|20i3 ). 
However, a significant increase in the young stellar popu- 
lations in the outskirts of the galaxies, for example, would 
radically change the FUV to optical SEDs, which are other- 
wise very well reproduced by our Maraston-mix (which does 
not include young stars). We are confident, therefore, that 
the SDSS spectra sample reasonably well the overall stel- 
lar populations in K+A galaxies. Barring significant stellar 
population gradients, the only free parameter in our simple 

5 some of the bhb models are tabulated for only two ages lOGyr 
and 15Gyr, making them much less convenient to use. 
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Figure 5. Left. SED typical of the K+A galaxies for which we have the complete FUV to MIR fluxes. The data are seen to be well 
approximated the MARASTON-mix models plus a Black-Body of T duat = 300°K. The blue triangles show the SPITZER 3.6/i, 4.5/i, 5.8//, 
8.0/i, and 24/t fluxes for this galaxy, which are seen to match within the errors the corresponding WISE fluxes. Right. Example of an 
SED that is not well fitted by the model at UV wavelengths. Both examples, however, show features that are typical of the whole sample: 
the Maraston models overestimate the near-IR fluxes; a single 300K blackbody is imposed on the Rayleigh-Jeans tail of the M05 models 
to account for the rising mid-IR flux . 



SED fitting approach is the zero point offset between the 
photometric and spectroscopic fluxes, which we calibrated 
out by matching the models to the SDSS r'-band fluxes. 



3.1 Emission-line Properties of K+A Galaxies 

We have seen that the overwhelming majority of K+A galax- 
ies are LINERs (i.e. "low ioniza tion nuclear em ission-line 
regions") originally proposed by iHeckmanl (jl98fj ) as a dis- 
tinct class of AGN. However, LINERs are actually a mixed 
bag of objects with very di fferent ionization mechanisms (e.g 
iTerlevich fc Melnicldll985r i . Thus, some LINERs are genuine 
AGN whereas most are actually powered by stars: they are 
'fake' AGN. 

CFS11 analyzed the spectra of a sample of several hun- 
dred thousand galaxies from the SDSS, more than half of 
which turned out to be LINERs. This is by far the most 
comprehensive study of the properties of LINERs we are 
aware of and contains an excellent up to date review of 
the 'LINER problem'. On the basis of extensive modelling, 
CFS11 showed that 'true' and 'fake' LINERs can be nicely 
separated using their WHAN diagram, and therefore pro- 
posed a new classification scheme for emission-line galaxies. 
This new scheme is presented in Figure that plots the 
WHAN diagram of K+A galaxies as in Figure [4] 

The new 'physically inspired' classification scheme of 
CFS11 divides emission- line galaxies into 4 groups as fol- 
lows: they separate LINERS into two groups, 'weak AGN' 
(wAGN), and 'retired galaxies' (RG). Thus, wAGN are ob- 
jects with 0.6 > EW(Hq) > 0.3A and log([NII]/Ha) > 



—0.4, which CFS11 claim to be powered by bona-fide 
nuclear black-holes. Conversely, RGs are objects with 
EW(Ha)<0.3A that are not powered by a nuclear engine, 
but by 'hot old low-mass evolved stars' (HOLMES) and 
therefore were actually 'fake' AGNs in the old classifica- 
tion scheme. As expected, the galaxies with strongest [Oil] 
& [OIII] lines lie closest to the AGN boundary as also do 
galaxies with 'significant' (i.e. >1.5%) young stellar popu- 
lations. In K+A galaxies the strong UV emission from the 
intermediate age populations may supply the UV flux which 
(in the HOLMES case) is actually provided by old Helium- 
burning stars, since the older stellar populations are not very 
massive in our galaxies. 

Although som e AGNs are known to exist in K+A galax- 
ies (|Liu et al.ll2007r ). Shin, Strauss & Tojeiro (2011) find that 
there are only a few bona fide AGNs associated with these 
galaxies and that the AGN is not related to the previous 
star formation episode. In our previous work, we showed 
that most K+A galaxies do not show evidence of AGN ac - 
tivity by stacking their 20cm signal l|Nielsen et al.l [20121 ). 
iBrown et al.l (|200fj| ') argue that AGNs are only responsible 
for quenchin g the more massive K+A galaxies, although 
ISnvder et all |2011i ) presents a model where K+A galaxies 
are the product of gas-rich mergers, where the K+A phase 
follows quenching by central AGNs (see also the discussion 
in Mendel et al. 2012). 

It is not obvious to us that most, or even some, of the 
objects in the wAGN zone are actually not 'fake' AGN. In 
order to probe further into this issue in the left panel of 
Figure [6] we show with different colors galaxies with partic- 
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Figure 6. Same WHAN diagram of Figure[4]but now using the new classification scheme proposed by CFS11. Thus, LINERs have been 
split into 'retired galaxies' (RG) and 'weak AGN' (wAGN), whereas star-forming galaxies (SF) now reside in the upper-left corner of the 
diagram. Strong AGN (sAGN) occupy the region previously ascribed to Seyferts. We have used different colors to code the strengths of 
[OII]3727 and [OIII]5007 in the left panel, and some relevant parameters from the STARLIGHT modes as indicated in the figure captions. 
wAGN are seen to have predominantly metal poor intermediate-age stars; lower masses; higher extinctions; and discernible young 
populations. 



ularly 'strong' (in this context!) [Oil] and [OIII] lines. These 
'strong' lined galaxies tend to concentrate in the wAGN 
zone, but many are also RGs, so we can't use these lines 
as criteria to distinguish true from fake AGNs. 

In the right-hand panel of Figure |S] we have color-coded 
the data according to some relevant parameters from the 
STARLIGHT models. The figure clearly shows that galaxies 
where the intermediate-age stars are predominantly metal- 
poor concentrate in the wAGN zone, while galaxies where 
the intermediate-age populations are metal rich tend to 
avoid that zone and concentrate in the region of RGs. Galax- 
ies with lower stellar masses (Mtot, higher internal extinc- 
tion, and the few galaxies with significant (i.e. >1%!) young 
stellar populations also tend to concentrate in the SF or 
wAGN domains. The above argues that the engine of the 
weak emission lines comes of the UV emission from inter- 
mediate age populations, rather than HOLMES as in other 
RGs studied by CFS11 

A priori, if wAGN are ionized by a central engine, we 
would not expect to find their emission-line properties to be 
correlated in any way with the stellar populations of their 
parent galaxies, in particular in the case of K+A galaxies 
where the populations of young stars are negligible. But we 
do see a correlation, which leads us to posit that the majority 
of the K+A galaxies in the wAGN zone may actually be fake 
AGN. 

The strong excess at A > 5[im observed in all K+A 
galaxies (e.g. Figures [5] and I12p imply that they have sub- 
stantial amounts of hot dust, although we do not know where 
this dust component is physically located within the galax- 



ies (see below for further discussion). The STARLIGHT ex- 
tinctions tell us that the dust appears to be well mixed with 
the gas, and that the intermediate-age population appears 
to be located in regions o f lower dust cont e nt. Th is is con- 
sistent with the model of IPoggianti fe~Wul |§000) where A 
st ars migrate out of n atal dust clouds, and the simulations 
of ISnvder et"ai1 i|201ll ) where the intermediate age popula- 
tion is revealed as dust is cleared by the AGN episode. We 
will return to this issue in the following sections when we 
analyze the ensemble properties of the dust as told by the 
SEDs. 



3.2 Gas-phase metallicities 

Before we leave the emission-line spectra of K+A galaxies, 
it seems relevant to examine whether the metallicities of the 
nebular component are similar or different to those of the 
dominant stellar populations. 

There are several semi-empirical methods to estimate 
the metallicities of HII regions for cases (such as ours) where 
the physical conditions of the ionized gas, and in particular 
the electron temperatures (T e ), are not well known. The 
most popular of th ese methods is the so-called -R23 origi- 
nally introduced bv lPagel et al.| || 19791) and later developed 
by nu merous studies (see e.g. iPena. Peimbert. fc Peimbertl 
12012] for a recent review). This method, however, is not well 
suited for our objects because it is ill-determined for metal 
rich nebulae, and because it requires the [OII]A3727,29 fluxes 
(-R23 = ( [Oil] + [OIII] ) /H/3) , which lie in the noisiest parts of 
our K+A spectra. 
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together with an extensive review of the literature are pre- 
sented by PP04 and we refer to that paper for further de- 
tails. A distinct advantage of both methods is that they are 
independent of extinction, although in K+A galaxies, the 
H/3 nebular line is heavily contaminated by the stellar ab- 
sorption from the A-stars, making it difficult to detect, and 
rather uncertain when actually detected. We are therefore 
forced to use Ha instead, and hence to correct the observed 
fluxes for extinction. 

This issue has been extensively discussed by CFS11, 
who derived an empirical relation between the nebular ex- 



tinction, A? 



estimated from the Balmer decrements in 



12 + log(0/H) 



their sample of hundreds of thousands of objects, and the 
extinction of the stellar spectra Ay resulting from their 
STARLIGHT models. CFS11 remark, as we do here, that their 
empirical relation between A v eb and Ay (probably) overes- 
timates the nebular extinction, which provided them with a 
suitable upper limit to the parameter space they were inter- 
ested in exploring. 

The opposite is true in our case where we are in- 
terested in a conservative lower limit for 03N2, which in 
any case does not depend very strongly on extinction. We 
therefore used the Ay values from STARLIGHT (without at- 
tempting to estimate the nebular component) to correct 
the observed fluxes for extinction. Following CFSI1, we 
assumed Case B recombination for the Balmer decrement 
(Ha/H/3=2.8), thus, 

03N2 = log(2.8 x [OIH]/[NII]) + 0.12 x Ay 

for the CCM extinction law. Figure [7] plots in the top panel 
03N2 vs. N2 for those galaxies in our sample with S/N>5 
in the continuum around Ha, and EW(Ha) > 0.1 (chosen to 
minimize contamination by spurious features). The dashed 
lines shows the limits of N2 > —0.5 and 03N2 < corre- 
sponding to nebulae with oxygen abundances 12+log(0/H) 
larger than solar (from Figs. 1&2 of PP04). PP04 also pro- 
vide analytical fits to estimate abundances as, 



Figure 7. Gas-phase abundances of K+A galaxies. The top panel 
shows a scatter plot of the two indices used as proxies for nebular 
metallicitiy: N2=log([NII/Ha) and 03N2=log([OIII]/[NII]). The 
lower panel shows the resulting abundance distributions using the 
best empirical calibrations from the literature. The dashed blue 
line shows the solar abundance of oxygen 12+log(0/H)=8.66. 



Fortunately, other empirical methods using the more 
benign [OIII]A5007 and [NIIJA6584 li nes have been exten- 
sively studied in the literature (e.g. iPettini fc P agel 2004 
henceforth PP04 and references therein), which are much 
better suited for our purposes, in particular because they are 
calibrated on a 'larger' (but still n <10) sample of metal-rich 
nebulae. Their calibration of these two proxies for metallic- 



12 + log(0/H) = 8.90 + 0.57 x N2 

and, 

12 + log(0/H) = 8.73 - 0.32 x 03N2 

shown in the bottom panel of Figure [JJ Within the uncer- 
tainties discussed above, the overwhelming majority of our 
K+A galaxies are seen to have nebular oxygen abundances 
well above the solar value (12+log(0/H)=8.66; P04). 

Notwithstanding the known discrepancy between the 
gas and stellar abundances in th e Orion region (which is 
now n arrowed to less than 0.1 dex - ISimon-Diaz k, Stasinskal 
2011), we are confident that the ionized gas and the 
intermediate-age stars in K+A galaxies have similar abun- 
dance distributions. 
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Figure 8. UV color-color diagram for 188 K+A galaxies with 
both FUV and NUV data from GALEX, corrected for foreground 
extinction, plotted in black. The blue points illustrate the locus of 
elliptical galaxies in the Coma cluster, and the red lines delineate 
the area occupied by horizontal-branch (HB) stars in Galactic 
globular clusters. The red points show integrated GALEX pho- 
tometry for Galactic GCs. The big green star shows the average 
colors of low-redshift QSOs and the arrow shows the effect of one 
magnitude of extinction for the CCM reddening law. 



3.3 UV fluxes 

Rising far-UV flux in early-type galaxies has been known 
to exist since the 1970s with the first space UV satellites, 
and it is now commonly attributed to the presence of hot 
horizontal branch stars (see O'Connell 1999 for a review). 
Figure [8] presents the colour-colour diagram commonly used 
in the literature to capture the broad statistical proper- 
ties of the ultraviolet SEDs of stars and galaxies. The blue 
points illustrate the locus of elliptical galaxies in Coma 
<|Smith. Lucev. fc Cartel I2012T I. These data have not been 
corrected for foreground extinction, which for Av = lmag 
and the CCM extinction law is represented by the arrow 
(the foreground extinction to Coma is of course much less 
than one magnitude). 

The red lines delimit the area occupied by HB stars 
in Galactic globul ar clusters (GC). adapted from the IUE- 
based diagram of lO'Connelj (|l999T) using the GALEX ob - 
servations of individual HB stars by ISchiavon et alj (|2012T ) 
to transform from IUE to GALEX, and the standard (g'- 
V) color of blue stars. The red points s how the integrated 
colour s of a sample of Galactic GCs from lDalessandro et al.l 
l|2012l ). Finally, the green star shows the average colours 
from a sample of 1600 SDSS QSOs at z ^ 0.4 for which we 
downloaded the relevant data from the GALEX and SDSS 
archives. 

K+A galaxies are not only substantially bluer than pas- 
sive ellipticals, as expected, but they are also significantly 
bluer than HB stars in GCs, consistent with the vacuum 



Figure 9. SED of one of the few K+A galaxy in our sample that 
has bluer FUV colors than typical QSOs. This also the only one 
in our sample with a significant young stellar population (35% of 
the mass) and lies well in the SF region in the WHAN diagram 
of Figure [4] However, it's SED is remarkably similar to a typical 
QSO. 



UV f luxes emerging from ~ 1 Gyr old stars l|Kavirai et al.l 
2008). K+A's are however not as blue as QSO's (green 
star), which are more than one magnitude bluer than even 
the bluest galaxies in our sample. There is one galaxy, 
SDSSJ142257.727+225441.43, however, that has UV colors 
typical of QSOs. This is actually one of the few galaxies in 
our sample that falls in the SF region of the WHAN di- 
agram (right panel of Figure and the only one with a 
significant population of young stars (39.5%). This is also 
one of the few galaxies in the sample with oxygen abun- 
dances lower than solar (Figure [7J), fully consistent with its 
STARLIGHT populations that are dominated by young and 
intermediate-age m etal poor stars, with no trac es of an old 
stellar component. IVeron-Cettv fc Veronl (120061) classify i t 
as a QSO based on the data by Darling fc Wegnerl l|l994f ). 
but this object shows no broad emission lines and no X-ray 
detection or radio signal in the FIRST survey. Interestingly, 
its SED shown in Figure [9] is surprisingly simi l ar to that of 
average radio-quiet QSOs (e.g. lAndreani et al.1 l|2003l )). with 
a prominent mid-IR excess at 10 — 20/i and a power-law con- 
tinuum in the blue. 



3.4 TP-AGB Stars 



iMarastonl (|2005t ) has suggested that the inclusion of TP- 
AGB stars in modelling intermediate age populations may 
dramatically alter their colours, derived ages, and stellar 
masses. Testing this hypothesis has proved difficult, as these 
stars contribute most of their light in the inf rared that is 
not ea sily a ccessible at hi g h reds hifts. However. iTonini et al.l 
(2009) and ITonini et all l^Old ) claim that models includ- 
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ing TP-AGB stars produce a better match to the colours 
of distant galaxies, while yielding ages and stellar masses 
in better agreement with the predictions of theoretical ex- 
pectations. Milner, Rose & Cecil (2011) obtain a good fit 
to the (CMD-calibrated) ages of M32 and NGC5102 using 
indices derived from infrared spectra and models including 
TP-AGB stars. lEminian et all (|2008l ) suggest that the in- 
frared colours of star-forming gal axies in the SDSS are better 
modeled by the iMarastonl (120051) models in cluding the TP- 
AGB contribution. iMelbourne et all |2012l ) claim that TP- 
AGB stars may contribute as m uch as 17% o f the A"-band 
flux in nearby galaxies. However. iKriek et al.l |2010l 'l find no 
evidence that TP-AGB stars are important in their sample 
of high er redshift post-starburst galaxies, while lZibetti et all 
l|2012h searched for the spectroscopic indices of TP-AGB 
stars in 16 z~0.2 K+A galaxies and fin d no strong evi- 
dence for their presence. The models by IConrov fe GunrJ 
il2010l) are in bette r agre ement with the older models by 
iBruzual fc Chariot! (|2003r) and yield a lower contribution 
from such stars. 

Figure [10] displays some of the broad-band colors com- 
monly used in the literature to characterize the SEDs of 
galaxies. (Notice that several of these colors have been 
shifted as indicated in the figure legend to minimize over- 
laps.) For each color we also plot the average color pre- 
dicted by our best-fit Maraston-mix plus the 300K black- 
body (which as discussed above we impose ad hoc to re- 
produce the observed mid-IR fluxes) shown by the arrows, 
and the dispersion of these predictions plotted as the ±ler 
horizontal bars on top of each arrow. 

The M05 models are seen in the figure to reproduce 
remarkably well the FUV to optical colors of K+A galaxies. 
This demonstrates that the FUV colors of K+A galaxies are 
bluer than HB stars because the intermediate-age (A) stars 
are actually quite luminous in the UV. However, the models 
systematically predict bluer K 3 —g' and K 3 — 3.4/i colors than 
observed. Figure [HI shows that these discrepancies actually 
arise because the models overestimate the fluxes in the NIR 
bands. 

As discussed in Section [2] we are confident that our 
photometric sample is not affected by systematic aperture 
matching effects. Nevertheless, we have argued that the 
influence of TP-AGB stars in the near-IR fluxes of high 
redshift galaxies is a contentious issue, which needs to be 
treated with extra care. Therefore, in Figure [TTJ we singled 
out galaxies classified as extended sources in the 2MASS or 
UKIDSS catalogues, for which we can be sure that the pho- 
tometric apertures are ideally matched. The figure shows 
that extended objects are indistinguishable from the rest of 
the sample proving that the NIR excess of the models is not 
due to systematics in the photometry. 

Surprisingly, our simplistic MARAST0N-mix+BB300K 
models appear to reproduce remarkably well the mid-iR col- 
ors of the sample, indicating that a single-temperature dust 
component suffices to reproduce the MIR excess in K+A 
galaxies (but see below). 

Another way of visualizing the systematic properties of 
the K+A galaxy SEDs is shown in Figure [T2l that presents 
the averaged SED of K+A galaxies binned in quartiles of 
the fraction of intermediate-age stars (inter) labelled with 
different colors. An Ensemble Maraston-mix+BB300K, ob- 
tained by averaging the stellar populations of all the galaxies 
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Figure 10. Distribution of observed colors of K+A galaxies 
using two different scales. The scale to the left corresponds to the 
visible to 3Afi colors while the scale to the right corresponds to 
the colours involving the FUV or 22fi fluxes. The arrows show 
the average values for the MARAST0N-mix+300K BB fits, and the 
dispersion of the sample indicated by the ±1<t bars on top of each 
arrow. A finer binning was used for (g' — r') and (K s — 3Afi) to 
avoid compressing the scale of the plot too much. 



in the sample with redshifts 0.04 < z < 0.3, is seen to over- 
estimate the observed NIR fluxes but underestimate the 3.4 
and 4.6 fluxes. This indicates that our models may be too 
simplistic to fully reproduce the SEDs. 

In particular, the discrepancy in the MIR would im- 
ply that either the temperature of the hot dust component 
changes significantly from object-to-object, which would be 
consistent with the observed dispersion in the WISE colors 
shown in Fig llOl or that more than one dust component is re- 
quired to reproduce the observations. However, some of the 
trends shown by Figure [12] may actually be due to features 
of the Maraston models for intermediate age stars. These are 
illustrated in Figure [T3] that shows the models averaged for 
stars between 0.9 — 2Gyr used in our Maraston-mix, color- 
coded by metallicity. 

The figure shows that the predicted UV fluxes depend 
dramatically on metallicity such that at A ~ 0.15/i metal 
poor stars are up to six orders of magnitude more luminous 
than metal-rich stars. This would predict that galaxies with 
the highest fractions of (intermediate-age) metal-poor stars 
should be systematically over-luminous in the UV. Since the 
models also predict (not shown here) that the UV fluxes 
of old metal poor stars are much larger than those of old 
metal-rich stars, reaching up to ~ 6 orders of magnitude at 
A ~ 0.15/i, the UV excess of galaxies with low fractions of 
intermediate-age stars, represented by the red symbols in 
Fig[l2] could actually be the telltale of HOLMES. 

In the NIR (magnified in the figure for clarity) the mod- 
els also show a significant dependence of metallicity. In par- 
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Figure 11. Ratios of predicted to observed near-IR fluxes in 
the H and Ks bands of 2MASS and UKIDS. The red symbols 
represent galaxies for which 2MASS isophotal and/or UKIDSS 
Petrosian magnitudes are available. The complete overlap of red 
and black dots shows that the observed ratios are not affected 
by aperture matching effects. The line for x=y indicates that the 
models overestimate H more than K. 



ticular, for Z=0.001 and Z=0.007, which approximately cor- 
respond to our metal-poor and metal-rich bins, the metal- 
poor model is the most luminous in the NIR and shows the 
strongest molecular bands, while the metal-rich model is the 
faintest and shows almost no molecular bands. All the other 
models (including Z=0.0001) lie in between and are roughly 
similar to one another. 

In the MIR the models systematically underestimate 
the observations at 3.4^i and (the 12p point is used 

to anchor the BB300K component), but the dust tempera- 
ture cannot be raised too much without compromising the 
fit at 22/i. However, the models have a discontinuous step 
at 2.5/x, which is a numerical feature introduced to extrap- 
olate the stellar models into the Rayleigh- Jeans regime at 
MIR wavelengths. As one would expect, this numerical step 
depends on metallicity, although the Z=0.001 step is sig- 
nificantly larger than for Z=0.0001. Figure [13] shows that 
for stars in our metal-poor bin (Z ~ 0.001) the 2.5/x step 
is about 0.5mag, which corresponds almost exactly to the 
discrepancy between the Ensemble Maraston-mix+BB300K 
and the observations in the MIR bands. 

We must conclude, therefore, that although the mod- 
els reproduce surprisingly well the FUV-to-optical colors 
of K+A galaxies, they provide no evidence that TP-AGB 
stars contribute to the infrared light in these p ost-starburst 
galaxi es, in agreement with t he results of iKriek et al.l 
(|2010t) and IZibetti et ail (|2012f l. The SDSS spectroscopy 
and the broad-band SEDs strongly in dicate that the 
older Bruzual fc Char otj (2003) and newer IConrov fc Gunnl 
(|201oF SSP models are adequate and that an increased con- 



Figure 12. Composite SED for K+A galaxies with varying frac- 
tions of intermediate age populations (indicated in the figure leg- 
end) compared with the normalized Ensemble MARASTON-mix ob- 
tained by averaging the STARLIGHT populations of all the galaxies 
in the sample fitted to the average fluxes plus a 300K black-body 
(black squares). In order to optimize the averages we have re- 
stricted the redshift range to encompass the bulk of the sample 
avoiding the extremes. Notice that the plot is in units of absolute 
luminosity, so the bars actually indicate the range of luminosities 
of the galaxies in our sample. Notice that the SEDs do not depend 
significantly on the intermadiate age fraction. 



tribution from TP-AGB stars is not necessary. In turn, this 
may have important implications for the nature of high 
redshift galaxies, whose colours and luminosities may then 
truly reflect thei r large ages, stellar ma sses and metallicities 
|Yan et al.ll2004l ; ICirasuolo et al.ll2010T iFl 



3.5 FIR excess 

Figure [T5] illustrates the fact that all our galaxies show sig- 
nificant excess at MIR wavelengths (in particular 12/i and 
22/x) compared to the M05 models. We have thus added 
a Black-Body of T dust = 300° K (BB300K) to describe the 
MIR upturn in the spectra. The temperature of the BB com- 
ponent was defined as discussed in Section \2. 31 and the nor- 
malization fixed by matching the 12/i fluxes as shown in 
Figure HQ We note h ere that the cluster objects studied by 
iDressler etHl (|2009T l show no significant excess in the MIPS 
band at 24/^m, possibly pointing to a real difference between 
cluster and field K+A galaxies. On the other hand, 6 out of 
the 18 galaxies they studied in Abell 851 have detectable 



6 The new generation of models by Claudia Maraston (private 
communication) having a smaller 2.5/x step will probably provide 
a better fit in the MIR, but the full models are not yet available 
for us to check whether they will also provide a better fit in the 
NIR 
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Maraston models averaged over: 0.9-2Gyr 
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Figure 13. Maraston SSP (Kroupa IMF) models for intermediate 
age showing the effect of metallicity, which is dramatic in the UV. 
The insert magnifies the near IR fluxes between 1 and 2.5fi. 



emission at (rest-frame) 16/xm, which is consistent with the 
fraction of K+A galaxies we detect in the reddest WISE 
band. 

There are various possibilities for the source of the 
mid infrared emission: dusty AGNs, dust enshrouded mas- 
sive stars a nd clusters, or circumst e llar dust from TP- 
AGB stars (IKelson fc Holdenl l201fj ; IDonoso et all |2012| ; 
IChisari fc Kelson! I2012T) . The 'typical' 22/x-MIR bump ob- 
served in K+A galaxies is actually quite sim ilar to the aver- 
age SE P of QSOs, as for example those by lAndreani et alj 
(2003), although the FUV up-turn is significantly steeper in 
QSOs. Of course the QSO may be so heavily obscured as to 
depress the UV contribution altogether. 

To gain further insight into the source of dust heating 
in K+A galaxies, Figure [14] shows the MIR excess index 
(12p — 3.4/j) (which is also a measure of the dust temper- 
ature) plotted as a function of internal extinction, Ay,int 
given by STARLIGHT. Note that we also used the STARLIGHT 
Ay to de-redden the observed fluxes, although we also veri- 
fied that our main conclusions do not depend on extinction. 

The figure shows a significant correlation between the 
MIR excess and the extinction seen by the stellar popu- 
lations. The slope of the correlation is -1.58+0.006, so the 
trend is highly statistically significant. The blue points show 
the galaxies with EW(Ha) >3A, which corresponds (except 
for just a few SF objects) to the locus of wAGN galaxies. 
Thus, for a given Ay wAGN have more MIR excess than 
RG. This may be interpreted as evidence that wAGN are 
enshrouded in dust that is not seen by the dominant 
stellar populations, and which must be rather hot in or- 
der to radiate at 12/i. 

An important feature of Figure \T%\ is that on average 
(black squares) K+A galaxies show a sizable FUV excess 
relative to the models; The right panel of Figure [14] shows 



that the UV excess, measured by (2270 — r'), is actually 
correlated with the temperature of the hot dust component, 
mea sured by (12/x — 3Au) Q 

IDonoso et all (|2012h argue that excess flux at 12 /im 
may originate from young stars and propose the 4.6-12 
color in WISE data as a measure of obscured star for- 
mation. Rising far a nd mid-IR fluxes ar e observed in ob- 
scured massive stars l|Whelan et al.|[201ll ) so a plausible in- 
terpretation of our results is that K+A galaxies are not 
truly 'dead' but that star formation continues within a 
dusty starburst. Dusty obscured starb ursts have been pro - 
posed as the origin of K+A galaxies bv lShiova et al.l l|200ll ). 
while limits to on - going star formation have been se t by 
Smail et all dl999T>; iMiller fc Owenl <|200ll) ; lGotol (|2004T ) and 



Nielsen et al 



(2012). 



On the other hand, dust heated by AGB stars 
in circumstellar envelopes may also become quite hot 
and migh t explain a mid/ f ar in fr ared excess, a s sug - 
gested by IKelson fc Holdenl (|201Cj) : IDonoso et all (|2012T ) : 
[Ch isari fc Kelsonl (|2012h . Theoretical understanding of this 
phenomenon does not allow us to fit models but the predic- 
tions suggest that the emission would peak at the interme- 
diate ages typical of K+A galaxies. 



4 SUMMARY 



CFS11 (|Cid Fernandes et alJlioTl have commented upon 
the challenges of conveying the information from large sam- 
ples of galaxies in an efficient and comprehensible way. Our 
sample of 'only' 808 K+A galaxies pales compared with the 
mega-studies undertaken by the SEAGal team, but we feel 
we have faced similar challenges. We have attempted the 
usual martingales of binning the data, or using plots in tech- 
nicolor, to convey the ensemble properties of K+A galaxies, 
but we still feel that it would be useful for the reader to 
recapitulate here the most important results. 

In summarizing these results, we will make an ontolog- 
ical leap of faith and generalize the properties of the 808 
galaxies in our sample to K+A galaxies as a whole. 

(i) K+A galaxies are significantly bluer in the FUV than 
Horizontal-Branch stars in Globular clusters; The FUV to 
optical colors of K+A galaxies are well reproduced by the 
M05 models as mostly arising from the intermediate-age (A) 
component; 

(ii) The overwhelming majority of K+A galaxies have 
weak but clearly detectable emission lines. About 75% of 
K+A galaxies are 'Retired Galaxies' where the nebular com- 
ponent is ionized by either intermediate-age stars, or by hot, 
old, low-mass, evolved stars (HOLMES in the nomenclature 
of CFS11), or both. The remaining ~15% fall in the cat- 
egory of 'weak AGN', considered by CFS11 to be genuine 
AGN; 

(iii) These wAGN have systematically stronger emission 
lines ([Oil], [OIII], Ha, & [Nil]), higher stellar extinction, 
lower fractions of intermediate-age metal poor stars, higher 
fractions of old metal poor stars, and higher fractions of 
young stars; 

7 The correlation is actually stronger for (1520 — g') albeit for a 
much lower number of galaxies 
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Figure 14. Left. Relation between infrared excess and the internal extinction in the galaxies derived from the STARLIGHT models 
showing a clear trend of 12)i fluxes increasing with extinction. This indicates that the dust component causing the IR excess is related 
to the dominant stellar populations. The blue symbols show that for a given amount of extinction, wAGN have systematically larger 
MIR excesses Right. UV excess plotted as a function of the (12/^ — 3.4/x) color, which for a single dust component is a proxy of dust 
temperature. Bluer colors indicate larger UV excess and higher dust temperatures. 



(iv) Compared to the M05 models, K+A galaxies do not 
show the predicted near-IR excess contribution from TP- 
AGB stars. It is likely that inclusion of these objects, at 
least at the levels postulated by the models and required to 
account for the red colors and lar ge masses of high redshift 
galaxies jTonini et al.ll2009l . l201oh , is not warranted; 

(v) At least one hot-dust component is required to repro- 
duce the fluxes at A > 5 /im. The MIR excess correlates 
with the stellar extinction indicating that at least part of 
the dust is mixed with the overall stellar populations. At a 
fixed extinction, wAGN and SF have stronger MIR excess 
than RG's indicating hotter dust temperatures. 

(vi) On average K+A galaxies show a significant FUV ex- 
cess compared to the models. The temperature of the hot 
dust component (as given by the MIR excess) correlates with 
the FUV excess, reinforcing our deduction that K+A galax- 
ies host massive dust enshrouded starbursts, which could 
possibly be the gen esis of bulges as recently suggested by 
iMendel et all (|2012| j. 

Although in the course of this investigation we have 
learned a great deal about K+A galaxies, we have not pro- 
vided a definitive answer to the fundamental question of 
whether K+A galaxies are currently forming stars. We found 
that the signpost of on-going star formation in these objects 
is not their blue FUV colors, which are mostly due to the 
intermediate-age (A) stars, but their MIR excess. In partic- 
ular, the excess appears to be well correlated with the 
WHAN type of the galaxies. This is illustrated in Figure [15] 
that shows the WHAN plot color-coded by the (12/i — 3.4/x) 
and (1520 — g') colors. This figure underlines our claim that 
wAGNs have stronger FUV excesses and hotter dust than 
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Figure 15. WHAN diagram for all 787 galaxies emphasizing ob- 
jects according to their 12 — 3.4 fim and FUV-g' colors. 



RGs, but also, most importantly, that many RGs also share 
the same properties. 

This, and the other 'technicolor' figures that we have 
displayed throughout this paper show that the spectral and 
photometric properties of K+A galaxies are continuous and 



16 J. Melnick & R. De Propris 



not multi-modal, as the WHAN partitioning could lead one 
to believe. In fact, a comprehensive review of the literature 
revealed no significant evidence for non-thermal radio or 
X-ray emission from K+A galaxies. This would support a 
model for K+A galaxies as havi ng evolved from dus ty star- 
bursts in the fashion described bv lShiova et al.l (|200lT ) , which 
may be continuing in highly obscured regions. 
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